Search for: All records

We present a set of six general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin accretion disks with different target mass accretion rates around black holes with spins ranging from nonrotating to rapidly spinning. The simulations use theM₁closure scheme with 12 independent frequency (or energy) bins ranging logarithmically from 5 × 10⁻³keV to 5 × 10³keV. The multifrequency capability allows us to generate crude spectra and energy-dependent light curves directly from the simulations without a need for special postprocessing. While we generally find roughly thermal spectra with peaks around 1–4 keV, our high-spin cases showed harder-than-expected tails for the soft or thermally dominant state. This leads to radiative efficiencies that are up to five times higher than expected for a Novikov–Thorne disk at the same spin. We attribute these high efficiencies to the high-energy, coronal emission. These coronae mostly occupy the effectively optically thin regions near the inner edges of the disks and also cover or sandwich the inner ∼15GM/c²of the disks.

 

Abstract We present and analyze a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a nonrotating, stellar-mass black hole. The simulations are initialized using the Shakura–Sunyaev model with a mass accretion rate of M ̇ = 3 L Edd / c 2 (corresponding to L = 0.17 L Edd ). Our previous work demonstrated that such disks are thermally unstable when accretion is driven by an α -viscosity. In the present work, we test the hypothesis that strong magnetic fields can both drive accretion through magnetorotational instability and restore stability to such disks. We test four initial magnetic field configurations: (1) a zero-net-flux case with a single, radially extended set of magnetic field loops (dipole), (2) a zero-net-flux case with two radially extended sets of magnetic field loops of opposite polarity stacked vertically (quadrupole), (3) a zero-net-flux case with multiple radially concentric rings of alternating polarity (multiloop), and (4) a net-flux, vertical magnetic field configuration (vertical). In all cases, the fields are initially weak, with a gas-to-magnetic pressure ratio ≳100. Based on the results of these simulations, we find that the dipole and multiloop configurations remain thermally unstable like their α -viscosity counterpart, in our case collapsing vertically on the local thermal timescale and never fully recovering. The vertical case, on the other hand, stabilizes and remains so for the duration of our tests (many thermal timescales). The quadrupole case is intermediate, showing signs of both stability and instability. The key stabilizing factor is the ability of specific field configurations to build up and sustain strong, P mag ≳ 0.5 P tot , toroidal fields near the midplane of the disk. We discuss the reasons why certain configurations are able to do this effectively and others are not. We then compare our stable simulations to the standard Shakura–Sunyaev disk. 

Relativistic magnetized jets, such as those from AGN, GRBs, and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and nonthermal radiation. Here, we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron–positron pair plasma. Generalizing recent treatments by Alves et al. and Davelaar et al., we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases, a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.